Peculiar long-range superexchange in Cu 2 A 2 07 (A = P, As, V) 
as a key element of the microscopic magnetic model 
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A microscopic magnetic model for a-Cu2P207 is evaluated in a combined theoretical and ex- 
perimental study. Despite a dominant intradimer coupling Ji, sizable interdimer couplings enforce 
long-range magnetic ordering at TV = 27 K. The spin model for a-Cu2P207 is compared to the 
models of the isostructural /3-CU2V2O7 and a-Cu2As207 systems. As a surprise, coupled dimers in 
a-Cu2P207 and alternating chains in Q-CU2AS2O7 contrast with a honeycomb lattice in /3-CU2V2O7. 
We find that the qualitative difference in the coupling regime of these isostructural compounds is 
governed by the nature of AO4 side groups: d-elements (A = V) hybridize with nearby O atoms 
forming a Cu-O-A-O-Cu superexchange path, while for p-elements (A = P, As) the superexchange 
is realized via 0-0 edges of the tetrahedron. Implications for a broad range of systems are discussed. 

PACS numbers: 71.20.Ps, 75.10.Jm, 75.30.Et, 75.50.Ee 
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I. INTRODUCTION 



The search for new magnetic ground states (GS) is an 
arduous challenge for condensed matter physics. Since 
about two decades ago, the rapid development of com- 
putational techniques has spurred theoretical studies 
on various magnetic models. The recent observation 
of magnetic monopoles, 1 anticipated by a theoretical 
prediction^ 2 - is one of the most remarkable examples ev- 
idencing the excellent potential and reliability of sound 
theoretical studies. 

Still, it is often tricky to find the way from a theo- 
retical result to its experimental confirmation. In par- 
ticular, the applicability of a theoretical model to a real 
material should be addressed. For magnetic systems, in- 
elastic neutron scattering (INS) experiments, which yield 
momentum-resolved excitation spectra for a wide range 
of energies and temperatures, are typically a method 
of choice. However, even the powerful neutron tech- 
niques are insecure against the ambiguous interpretation 
of experimental results. JThus, for instance, the analy- 
sis of INS data in Ref. y correctly disclosed the irrel- 
evance of the two-leg Heisenberg spin-ladder model for 
the ambient-pressure modification of (VO)2P2C>7. At 
the same time, the experimental spectra exhibited a pro- 
nounced high-energy mode, and the authors of Ref. [j in- 
terpreted it as a signature of a two-magnon bound state. 
Subsequent studies based on alternative techniques re- 
futed this conjecture and revealed that the high-energy 
mode and the double-gap feature arise from two inequiva- 
lent alternating Heisenberg chains (AHC) present in this 
system^ As another example, the frustrated Heisenberg 
chain LiCu 2 02, had been proposed to imply two relevant 
antiferromagnetic (AF) exchange integrals J\ and J2, 
based on thermodynamical measurements and neutron 
diffraction data,— while the combined nuclear magnetic 
resonance (NMR) and density functional theory (DFT) 
study (Ref. @) as well as the INS investigation (Ref. 0) 
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FIG. 1. (Color online) Crystal structure and exchange cou- 
pling paths in CU2A2O7 (A = P, As, V). Squares and triangles 
represent CUO4 plaquettes and AO4 tetrahedra, respectively. 



revealed the ferromagnetic (FM) nature of J\, 

Recently, density functional theory (DFT) calculations 
developed into an appealing alternative, complementary 
to experimental techniques. At present, the performance 
and reliability of DFT codes allow us to evaluate quanti- 
tative magnetic models even for complex materials, hav- 
ing only crystallographic information at hand. More- 
over, DFT calculations are capable of even improving 
the structural input,— thus overcoming the problem of 
weakly scattering elements (for instance, H for X-rays or 
V for neutrons) and the subsequently large error bars in 
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the experimental atomic positions. 

The inverse problem — to find a crystal structure for 
a specific spin lattice — is even more complex. Here, 
semiempirical Goodenough-Kanamori- Anderson (GKA) 
rules, 9 formulated about GO years ago, still provide in- 
valuable support. However, these rules are at best not 
universal. In particular, the GKA rules cannot be ap- 
plied to long-range couplings that involve more than one 
ligand atoms. The abundance of such couplings in real 
materials strongly motivates the development of further 
rules to account for these situations. 

Semiempirical rules are especially helpful for a fine tun- 
ing of the material, because chemical substitutions usu- 
ally alter exchange couplings or even the spin lattice. For 
example, the substitution of Br for CI changes the singlet 
ground state of (CuCl)LaNb 2 07 toward columnar AF or- 
dering in the bromine compound.— The replacement of 
Zn by Mg in kapellasite Cu3Zn(OH)gCl2 leads to hay- 
deeite Cu3Mg(OH)eCl2 with a reduced nearest-neighbor 
(NN) coupling and modified spin physics. 11 

In the present study, we perform a comparative investi- 
gation of the isostructural CU2A2O7 (A = P, As, V) com- 
pounds and consider an instructive example that is rele- 
vant for a family of systems with tetrahedral anions. The 
pentavalent nature of the A atoms ensures the magnetic 
3d 9 electronic configuration of Cu 2+ . The thortveitite- 
type crystal structure of CU2A2O7 (Fig.[IJ top) contains 
magnetic layers of Cu 2 06 structural dimers (two edge- 
sharing CUO4 plaquettes) coupled by non-magnetic AO4 
tetrahedra. The neighboring layers are connected by two 
AO4 tetrahedra forming A2O7 pyrophosphate, pyroarse- 
nate or pyrovanadate groups (Fig. [TJ top right). 

In this study, we disclose the underlying reason for 
a variety of magnetic behaviors — coupled spin dimers 
in C1-CU2P2O7, alternating spin chains in a-Cu2As207, 
and a honeycomb spin lattice in /3-CU2V2O7 — realized 
within the same basic crystal structure. Using DFT cal- 
culations, we evaluate five relevant exchange couplings 
(Fig. [D bottom). Surprisingly, these couplings are dras- 
tically affected by the electronic state of the side groups' 
central atom A, whether it is a p (P, As) or a d (V) el- 
ement. Understanding the crucial role of non-magnetic 
AO4 groups is the key for directed chemical substitution, 
which may open the door for a tunable magnetism. 

The magnetic properties of /3-CU2V2O7 have been 
actively studied*^ Very recently, we derived the ac- 
tual spin model and identified an almost ideal spin-1/2 
honeycomb lattice in this compound. 13 Considering the 
crystal structure of /3-Cu 2 V 2 07, this magnetic model 
looks highly non-trivial. Puzzled by the driving force of 
such an unusual magnetic coupling regime, we prepared 
the isostructural system a-Cu2P207 and investigated its 
magnetic behavior to allow for a comparative study. 

CU2P2O7 has two polymorphic modifications. The 
low-temperature (a) and the high-temperature (/?) 
phases reveal similar monoclinic crystal structures 
(Fig.fTJ), yet the j3- modification shows a twice smaller 
c lattice parameter and disordered arrangement of the 



oxygen atom shared by the two PO4 tetrahedra of the 
pyrophosphate group*^ In the following, we consider 
the a-modification which is fully ordered and stable be- 
low 350 K. We would like to emphasize that this a- 
modification is isostructural to the /3-modification of 
CU2V2O7. The magnetic behavior of a-Cu2P207 was 
studied in the early 1970s. An AF ordering was found 
at 27 K, with tentative models of the magnetic structure 
established from neutron diffraction^ and NMRJ^ How- 
ever, no information on the microscopic magnetic model 
has been reported. 

This paper is organized as follows. We introduce the 
methods used in our study in Sec.|H] In Sec. lIIH we 
present our experimental as well as theoretical investiga- 
tion of a-Cu 2 P207. In Sec. lIVl we compare the magnetic 
models of Cu 2 ^4207 (A = P, As, V), and discuss the origin 
of their essentially different nature. Finally, in Sec. El we 
summarize the results and give a short outlook. 



II. METHODS 

A single-phase powder sample of a-Cu2P 2 07 was pre- 
pared by a solid-state reaction of CuO and NH4H 2 P04. 
A stoichiometric mixture of the reagents was fired at 
250 °C for 6h in air to remove ammonia and water. 
The sample was further placed into a sealed quartz tube 
and heated at 800 °C for 12 h. The resulting light-gray 
powder was analyzed by X-ray diffraction (Huber G670 
Guinier camera, CuK a i radiation, image plate detec- 
tor, 3-100° angle range). A Rietveld refinement con- 
firmed the formation of a-Cu2P207 and did not show 
any traces of impurity phases (in particular, the ad- 
mixture of /?-Cu2P 2 07 can be ruled out). In contrast, 
a high-temperature annealing in air always produced a 
Cu3(P04)2 impurity. 

Magnetic susceptibility (x) was measured with a 
SQUID magnetometer (Quantum Design MPMS) in the 
temperature range 2-380 K in applied fields up to 5T. 
High-field magnetization measurements (at T=1.4K) 
were performed at the Dresden High Magnetic Field Lab- 
oratory in pulsed fields up to 60 T. Details of the experi- 
mental procedure are described elsewhere^ 

The electron spin resonance (ESR) spectra of a pow- 
der sample of a-Cu2P207 were recorded at X-band fre- 
quencies (9.4 GHz) for temperatures 5-300 K. The spec- 
tra were fitted with a powder average of two Lorentzian 
lines corresponding to two (effective) components of the 
g-tensor in axial symmetry: g± and g>\ . 

DFT calculations were performed using the full- 
potential code fplo9 . 00-31 M Within the local (spin) 
density approximation (L(S)DA), the exchange and cor- 
relation potential of Perdew and Wang has been usedJ^ 

For the LDA calculations, we chose a mesh of 1008 
k points (207 in the irreducible wedge). Spin-polarized 
supercell LSDA+C7 calculations were performed on 12, 
18, and 32 k points meshes^S All calculations were done 
for the experimental crystal structure from Ref. [T3 . To 
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address the consistency of the structural data, we calcu- 
lated forces and cross-checked the LDA results by using 
the generalized gradient approximation (GGA). Since the 
resulting forces do not exceed leV/A [the mean value is 
0.7(3) eV/A in both LDA and GGA], no further relax- 
ation has been performed. The calculations were care- 
fully checked for convergence. 

Quantum Monte-Carlo (QMC) simulations were per- 
formed on 24x24 finite lattices with periodic boundary 
conditions using the looper 2 * and dirloop_sse (directed 
loop in the stochastic series expansion representation) 22 
algorithms from the ALPS package.— 



III. MAGNETISM OF a-CU 2 P 2 7 
A. Magnetization and ESR measurements 

The magnetic susceptibility x(^) (Fig- [U top) ex- 
hibits a broad maximum at T max = 54 K indicative of 
low-dimensional behavior. The Curie-Weiss (CW) fit 
yields an AF Weiss temperature 0cw = 61K and the ef- 
fective moment of 1.89 /ie, slightl y larger t han the spin- 
only effective moment /i c g = 2y/S(S + 1) « 1.73 [Xb for 
5 = 1/2. The deviation from the CW fit below - 140 K 
originates from the onset of AF correlations that lead to 
a long-range magnetic ordering at low temperatures. 

The magnetic ordering transition is evidenced by a 
weak kink in the low-field x(T) curves around Tjv — 27 K 
(circles and solid line in the upper inset of Fig. [2|). The 
data measured at 5 T reveal a well-defined kink at the 
same temperature. The difference between the data col- 
lected in low and high fields below Tjv should be related 
to the spin-flop transition around 1 T. The susceptibil- 
ity upturn below 10 K is likely caused by a paramagnetic 
impurity contribution. 

In order to estimate the <?-value and prove the AF or- 
dering in CU2P2O7, we measured the electron spin reso- 
nance (ESR) of the Cu 2+ spins as a local probe (Fig [3]). 
The temperature dependencies of the integrated ESR ab- 
sorption and x(T) agree with each other (not shown). At 
elevated temperatures, both components of the <?-tensor 
(g± and gn) show typical values for the Cu 2+ ion. An 
abrupt change in g± and g>\ around T/v = 27K indicates 
the onset of a strong internal magnetic field and thus 
confirms the magnetic ordering transition observed in the 
X(T) data. 

The magnetization M(H) curve (Fig. [2j bottom) was 
scaled using the low-field (MPMS) data. It lacks any 
signatures of a spin gap and exhibits a linear behavior, 
incompatible with a simple dimer magnetism. 
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FIG. 2. (Color online) Top: magnetic susceptibility x(T) °f 
CU2P2O7 together with the Curie- Weiss (CW) and the quan- 
tum Monte-Carlo (QMC) fits (see text). Inset: field depen- 
dence of x(T) Bottom: high-field magnetization of CU2P2O7 
with a QMC fit. Inset: scaling of the high-field magnetization 
curve with the MPMS data. 
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B. DFT calculations 

To evaluate the leading exchange couplings, we per- 
form DFT band structure calculations. For a-Cu2P207, 
LDA yields a valence band with a width of 10 eV (Fig. |H 



FIG. 3. (Color online) Temperature dependence of the g- 
factor measured by the ESR of Cu 2+ . The inset shows a typ- 
ical derivative ESR spectrum (circles) fitted by two powder- 
averaged Lorentzian shapes (solid line). 
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FIG. 4. (Color online) Density of states (top), band structure 
(bottom, left), and orbital-resolved density of states (bottom, 
right). The bands with 3d x 2 _ y i character are shown bold 
(green). The Fermi level is at zero energy. 



top), slightly larger than typical values for cuprates (7- 
8eV)fiii2i The density of states (DOS) exhibits a peak 
at the Fermi level £f, evidencing a metallic GS, in con- 
trast to the expected insulating behavior. This drawback 
of the LDA originates from the underestimation of cor- 
relation effects, intrinsic for the 3d 9 electronic configu- 
ration of Cu 2+ . To account for the correlation effects, 
we use two complementary approaches: (i) mapping the 
LDA bands onto a Hubbard model and (ii) adding cor- 
relations in a mean-field way within the band structure 
scheme (LSDA+f/ calculations). 

Despite the problem of correlations, LDA yields reli- 
able information on the relevant orbitals and couplings. 
First, we project the DOS onto local orbitals of a Cu04 
plaquette. We find that the Cu states at ef have 3d x 2 _ y i 
character (Fig.Sl bottom right), resolving it as the mag- 
netically active orbital^ Further on, the relevant cou- 
plings were identified using Wannier functions (WF) for 
the Cu 3d x 2_ y 2 states. 26 The resulting fit is in excellent 
agreement with the LDA bands (Fig. 21 bottom left) and 
yields five relevant couplings. Two of them — the cou- 
plings within the structural dimer t\ and between the 
structural dimers t\ (Fig. [TJ middle) — form alternating 
chains along the b axis. The chains are connected by two 
types of interchain couplings: t- lc i and i; C 3 run through 
a single AO4 tetrahedron, whereas t- 1C 2 connects two Cu 
atoms with a double bridge of AO4 tetrahedra (Fig. [TJ 
bottom). Together with t\ and t' Xl these interchain cou- 
plings form the magnetic layers, parallel to (101), as 
shown in Fig. [T] The interlayer coupling of 30meV is 
realized via pyrophosphate groups P2O7. Since it is sub- 
stantially smaller than the leading terms, we neglect this 
coupling in the minimal model. 

Here, we briefly remind the reader of the essential re- 
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TABLE I. Transfer (ti, in meV) and exchange (Ji, in K) in- 
tegrals for Q-CU2P2O7 and /3-CU2V2O7. J AF (in K) are cal- 
culated as Atf /UcB (see text). 
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suits for /3-CU2V2O7 in Ref. HI the WF analysis of the 
LDA bands revealed that t\ , t[ , and t ic2 couplings are rel- 
evant, while other terms are negligibly small (Table IJ). 
A mere comparison of ti for ev-Cu2P207 with the respec- 
tive values for /3-CU2V2O7 readily yields several impor- 
tant results. First, the coupling t\ within the structural 
dimer is essentially the same, while the interdimer cou- 
pling t[ is slightly enhanced in CV-CU2P2O7 compared to 
/3-CU2V2O7. The similarity of the intra-chain terms is in 
sharp contrast to the profound change in the interchain 
coupling regime: a-Cu2P207 has sizable £; c i and i; C 3 with 
the negligibly small ti c2 , while for /3-CU2V2O7 the sit- 
uation is exactly the opposite. This results in clearly 
different band dispersions in a-Cu2P207 (Fig. [4]) and /3- 
CU2V2O7 (Fig. 3 in Ref. [H). 

To obtain the values of AF exchange integrals J AF , we 
map the transfer integrals ti from the WF-based anal- 
ysis onto a Hubbard model. For the half-filling and 
the strongly correlated limit U c ff^>ti, both well justi- 
fied for undoped cuprates, magnetic excitations can be 
described by a Heisenberg model, while J 2 AF is expressed 
as J AF =4t 2 /[/ eff . This way, using a realistic value of 4.5 
eV for the Coulomb repulsion J7 e ffr^ we obtain the AF 
part J AF of the total exchange Ji (Table Q}. 27 

For the NN exchange coupling path Xi, the Cu-O- 
Cu angle is close to 90° (98.7° for /3-Cu 2 V 2 7 , 100.4° 
for a-Cu2P207), suggesting a sizable ferromagnetic con- 
tribution according to the GKA rules £ Indeed, in many 
related systems this leads to a sizable reduction in the 
respective exchange Ji^- Moreover, in case the FM con- 
tribution exceeds </ AF , the resulting Ji exchange becomes 
FMiS Therefore, to complete the spin model, it is neces- 
sary to calculate the values of total exchange integrals 
Ji, being the sum of the AF and FM contributions. 
We evaluate the numerical values by performing spin- 
polarized LSDA+C/ supercell calculations for various 
collinear spin arrangements and map the total energies 
onto a classical Heisenberg model. Following the conclu- 
sions of Ref. [TI], we choose the around-mean-field double- 
counting-correction (DCC) and Ud = 6.5 eV^°- The re- 
sulting exchange integrals Ji are given in Table 

For a-Cu2P207, a sizable FM contribution of jf M ~ 
—230 K considerably reduces the coupling within the 
structural dimer, as in /3-CU2V2O7. Again similar 
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Cu 2 P 2 7 C112AS2O7 Cu 2 V 2 7 






FIG. 5. (Color online) Top: Wannier functions for Cu 3d x 2 _ y i states with relevant interchain couplings (lines). Bottom: spin 
models for Cu2^420 7 compounds. The thickness of a line reflects the magnitude of the respective coupling. Filled and empty 
circles represent up and down spins in the antiferromagnetically ordered state. 



to /3-CU2V2O7, the interdimer coupling J[ has a tiny 
FM contribution only. 3 - However, in a-Cu2P207, J[ 
is clearly dominant and seemingly favors a dimcr-likc 
magnetism with magnetic dimers located between the 
structural dimers Cu2 06- All interchain couplings are 
long-range, and their FM contributions are likely to be 
small. Indeed, we find J™«J™~-20K consistent 
with J™ = -10 K found for /3-Cu2V 2 7 . 13 Noteworthy, 
for the closely related system Cu2(P03)2(CH2), simi- 
lar values have been reported.— To summarize, for a- 
CU2P2O7 we arrive at the two-dimensional (2D) mag- 
netic model with a strongly dominant J[ and three ad- 
ditional relevant couplings, rather different from the 2D 
honeycomb lattice system /3-CU2V2O7 as well as from 
the structurally-related but essentially one-dimensional 
(gapped Ji — Jj chains) (VO)2P20 7 in its high-pressure 
modification.— 

Despite the 2D nature of the spin models, the mag- 
netic GSs of a-Cu2P 2 7 and /3-Cu 2 V 2 7 are AF- 
orderedi 15 ' 35 Since thermal fluctuations in purely 2D 
Heisenberg systems exclude long-range magnetic order at 
finite temperatures,— the coupling perpendicular to the 
magnetic layers should be addressed. In /3-Cu2V20 7 , ef- 
fective superexchange paths along the [V2O7] pyrovana- 
date groups give rise to a sizable long-range interlayer 
exchange J± = 17KJ£ In contrast, the interlayer cou- 
plings in CV-CU2P2O7 are much weaker and short-range: 
although the transfer integral associated with the short- 
est interlayer coupling (Cu-Cu distance 3.33 A) amounts 
to t±i =31meV, the LSDA+[7 calculations reveal a siz- 
able FM contribution that largely reduces this coupling 



to about J±i = 1 K. The next largest (t±2 — 18 meV) in- 
terlayer coupling corresponds to the Cu-Cu distance of 
3.88 A, and gives rise to a magnetic exchange of 3-4 K. 

Using the interlayer couplings resulting from the DFT 
calculations, we find one more difference between a- 
CU2P2O7 and /3-CU2V2O7. In /3-CU2V2O7, the magnetic 
structure should break the C-centering of the unit cell— 3 - 
whereas the long-range order in a-Cu2P20 7 inherits the 
full crystallographic symmetry (the C2/c space group). 
The distinct patterns of the magnetic ordering should be 
resolvable by high-precision neutron diffraction and could 
be an ultimate experimental test of our predictions. The 
lack of magnetic reflections in the available neutron data 
for a-Cu2P207 (Ref. ITsT) indicates that the magnetic or- 
der does not break any crystallographic symmetry, nor 
leads to the formation of a magnetic superstructure, in 
accord with our microscopic model. 



C. Model simulations 

To justify the proposed spin model for a-Cu2P207, we 
perform QMC simulations of the magnetic susceptibil- 
ity and the high-field magnetization data. First, to fit 
the calculated susceptibility to the experiment (Fig. [2]) , 
we keep the ratios Ji/J(=0.33, J; Cl / J[ = 0.45, and 
J ic JJ[ = 0.40 from the DFT calculations (Table [J), and 
vary the absolute values of J[ and gM This way, we ob- 
tain an excellent agreement with the experimental curve 
down to the ordering temperature T^ — 27K (Fig. 
The resulting J[ — 79K and g = 2.21 are in reasonable 
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agreement with J[ = 102 K from the LSDA+J7 calcula- 
tions and the powder-averaged g={2g±_ + g , ii)/3 = 2.18 
from ESR. 

The simulated magnetization curve (Fig. [5J bottom) is 
scaled adopting J[ — 79 K and g = 2.21 from the fit to the 
magnetic susceptibility. The overall agreement with the 
experiment is rather good, while the discrepancy likely 
arises from the imperfectness of the scaling of the experi- 
mental high-field curve. In particular, the scaling is done 
based on a very narrow low-field region (up to 5 T, see the 
inset of Fig. [2]) , which might be affected by paramagnetic 
impurities^ 



IV. DISCUSSION 

The spin models of a-Cu2P20y and /3-CU2V2O7 are 
essentially different: whereas a-Cu2P20y is a system of 
coupled dimers, /3-CU2V2O7 is an almost isotropic hon- 
eycomb lattice. In /3-CU2V2O7, a strong FM contribu- 
tion practically "switches off" the coupling within the 
structural dimers (Ji), leaving only the couplings be- 
tween the structural dimers ( J[) and between the chains 
(Ji C2 ) relevant (Fig.O bottom right) On the contrary, 
a-Cu2P207 has a dominant J{, smaller J lCl and J lC3 (see 
Fig. [U for the notation of the exchange couplings) as well 
as even smaller Ji (Fig. [SI bottom left). The striking dis- 
similarity of the two models leads to a question, whence 
this difference originates. 

To address this puzzling issue, we consider another 
Cu2^207 system, a-Cu2As207, and briefly report the 
key results of our DFT calculations. In particular, the 
WF analysis yields values of transfer integrals f j similar 
to those of a-Cu2P207. The total exchange integrals, 
derived from the results of LSDA+C7 calculations, are 
somewhat larger than the respective exchange couplings 
for a-Cu2P207. However, the only pronounced differ- 
ence between the magnetic models of a-Cu2P207 and a- 
CU2AS2O7 is the strong enhancement of the NN coupling 
J 1 in the latter system (Fig. [5J bottom middle). A de- 
tailed experimental and theoretical analysis of the mag- 
netic properties of a-Cu2As207 is presently underway^ 

So far, we find the resemblance of the spin models of 
ct-Cu2P207 and a-Cu2As 2 07, whereas /3-CU2V2O7 is re- 
markably different. The common feature of these com- 
pounds is the sizable AF coupling J[ forming spin dimers. 
The intcrdimer couplings are, however, different. In a- 
CU2P2O7 and a-Cu2As207, the interdimer couplings run 
through a single ^404 tetrahedron (J; Cl and Ji C3 ), while 
in /3-CU2V2O7 only the coupling via a double bridge of 
AO4 tetrahedra is relevant. 

As we know from the analysis of the band structure, 
the magnetism of Cu2^207 systems originates from the 
half- filled Cu 3d x 2 _ y 2 orbitals. Above, we used the WF 
for this magnetically active orbital to obtain the relevant 
couplings (Sec. IIII B|) . However, the fact that WF are 
defined in real space, makes them capable of tracing the 
superexchange path4S In Fig. [5J we plot the WF for the 



three systems. 

As expected from the similarity of the band struc- 
tures, the WFs for a-Cu 2 P207 and a-Cu2As207 are very 
similar. 41 Note (i) the absence of P or As contributions 
and (ii) the antibonding combination formed by WF con- 
tributions along Cu-O-O-Cu paths, favoring t\ Q \ and £i C 3 
couplings. In contrast, the WF of /3-CU2V2O7 shows siz- 
able weight on V atoms, exhibiting a shape reminiscent 
of the d 3z 2_ r 2 orbital. These V states hybridize with 
the nearby O orbitals, thus forming an effective superex- 
change path Cu-O-V-O-Cu. The V-0 hybridization 
drastically alters the interchain coupling regime, now fa- 
voring J- lC2 exchange via a double bridge of VO4 tetrahe- 
dra. 

As evidenced by numerous examples, the replacement 
of anionic groups can stabilize different crystal struc- 
tures. Thus, the replacement of Se by a larger Te 
atom transforms the chain- like structure of CuSe20s 
into complex 2D layers of CuTe2C>5. Not surprisingly, 
this structural change drastically affects the magnetism: 
while CuSe2C>5 shows Heisenberg chain physics and AF 
order^ CuTe205 has a dimerized magnetic GS4S How- 
ever, it is commonly believed that the replacement 
of a nonmagnetic part, not accompanied by structural 
changes, should have a minor impact on magnetism. Our 
study shows that this conjecture is not universal. In par- 
ticular, the electronic configuration of the nonmagnetic 
group should be addressed. Thus, the d-states of V in j3- 
CU2V2O7 are stron gly mixed with the p states of oxygen 
(see Fig. 2 in Ref. Il31 ). allowing for a strong superex- 
change along the largely covalent V-0 bonds. On the 
contrary, p-cations like P or As donate electrons to the 
neighboring O atoms, giving rise to more ionic P(As)-0 
bonds. In this case, electron hopping between P(As) and 
O atoms becomes energetically unfavorable, thus confin- 
ing the superexchange to the Cu-O-O-Cu path. 

Regarding the superexchange in transition-metal phos- 
phates, the lack of the P contribution to the magnetic 
orbital and the crucial role of M-O-O-M pathways has 
been recently established for the related spin- -i V +4 com- 
pounds. A systematic study of model structures with dif- 
ferent bond distances and angles puts forward the M-0 
and 0-0 distances as key structural parameters that de- 
termine the coupling strength . 43 ! 44 Another relevant fea- 
ture is the number of bridging tetrahedral groups. Given 
similar bond distances, a double bridge makes two su- 
perexchange pathways, which are more efficient than a 
single bridged This trend holds for a-Cu 2 P207, where 
the double-bridge coupling J[ is stronger than the single- 
bridge exchanges J- lc i and J- 1C 3 (Table [TJ> . An important 
difference between copper and vanadium phosphates is 
the stronger metal-oxygen hybridization in the former 
case. The hybridization, amplified by the cr-type bonding 
between Cu d x 2_ y 2 and oxygen p orbitals, leads to a non- 
negligible contribution of second-neighbor oxygen atoms 
to WF (see Fig. ©. In BiCu 2 P0 6 , such "tails" of WF 
are responsible for the pronounced difference between the 
seemingly equivalent long-range couplings^ A compar- 
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ative study of Cu +2 and V +4 compounds could facilitate 
further development of empirical rules for the superex- 
change. The CU2A2O7 structures are especially apposite 
for this purpose due to their similarity to (VO)2P207, 
CsV 2 5 , VOSe0 3 , and other V+ 4 oxides. 

Generalising the results, our study evidences that for 
the tetrahedral non-magnetic ^40*4 groups with the cen- 
tral p-atom the magnetic superexchange is realized via an 
0-0 edge of the tetrahedron. In contrast, a tetrahedrally 
coordinated 3d element gives rise to a strong hybridiza- 
tion with O 2p orbitals and a concomitant change in the 
superexchange path. This finding constitutes an empiri- 
cal rule to estimate the relevant long-range couplings in 
a wide range of low-dimensional systems containing non- 
magnetic tetrahedral side groups. 



V. SUMMARY AND OUTLOOK 

To summarize, we presented a comprehensive exper- 
imental (magnetic susceptibility, ESR, and high-field 
magnetization) and theoretical (full-potential DFT cal- 
culations) investigation of the magnetism of a-Cu2P20y. 
This material implies a two-dimensional magnetic model 
of strongly coupled magnetic dimers, which interestingly 
does not lead to a singlet ground state: sizable interdimer 
couplings enforce the antiferromagnetic ordering at 27 K. 

The coupled dimer model of a-Cu2P207 is essentially 
different from the honeycomb lattice model, realized 
in the isostructural system /3-CU2V2O7. Based on a 



comparison between a-Cu2P207, a-Cu2As207, and /?- 
CU2V2O7, we unravel the crucial importance of the non- 
magnetic AO4 side groups for the magnetism of Cu2^207 
(A = P, As, V) systems: whereas a d-atom (V) strongly 
hybridizes with the ligand O orbitals, the p-states (P, As) 
do not contribute to the superexchange path. Within the 
structure of Cu2^42 07 compounds, this leads to a drastic 
change in the topology of the interchain couplings. The 
forthcoming experimental study of a-Cu2As207 should 
provide a deeper insight into the magnetism of this sys- 
tem and challenge our findings. Our results can be ap- 
plied to tune the magnetism of systems comprising a 
similar structural motive: for instance, Cu2(P03)2(CH 2 ) 
(Ref.Hl or (VO) 2 P 2 07 (Ref.|I|. 

Our study emphasizes the danger to use isostructural 
compounds as reference or analog. Such comparisons, 
unless based on a precise knowledge of the underlying 
electronic structure, are at best not justified and can be 
completely misleading. Instead, they are amenable to a 
careful verification. The latter can be reliably performed 
by electronic structure calculations that provide a micro- 
scopic insight into the electronic and magnetic properties 
even of complex compounds. 
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